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Abstract: A path within the parameter space of detuning and pump 

power is demonstrated in order to obtain a single cavity soliton (CS) with 
certainty in SiN microring resonators in the anomalous dispersion regime. 
Once the single CS state is reached, it is possible to continue a path to 
compress it, broadening the corresponding single free spectral range (FSR) 
Kerr frequency comb. The first step to achieve this goal is to identify the 
stable regions in the parameter space via numerical simulations of the 
Lugiato-Lefever equation (LLE). Later, using this identification, we define 
a path from the stable modulation instability (SMI) region to the stable 
cavity solitons (SCS) region avoiding the chaotic and unstable regions. 
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1. Introduction 


An optical frequency comb is a light source with a number of highly resolved and nearly 
equidistant spectral lines. Since its introduction, multiple important applications have been 
demonstrated in areas such as communications, metrology, spectroscopy, astronomy and 
optical clocks. Optical frequency combs can be generated using mode-locked lasers or 
electro-optic modulation of continuous-wave light. Since 2007, multiple experiments have 
reported optical frequency comb generation by means of Kerr nonlinearity (wave mixing 
process) in microresonators, which offer potential for highly compact and portable solutions. 
Such combs are termed Kerr optical frequency combs or simply Kerr combs CHID- 

The understanding of underlying processes and dynamics in Kerr comb generation is 
critical to move this technology further to industry and applications. Chembo and Yu in JT2j 
introduced one of the first simulation approaches using modal expansion. More recently, the 
Lugiato-Lefever equation (LLE) 1131 has been widely adopted 114-24]|. Furthermore, Chembo 
and Menyuk in [15J demonstrated the equivalence between the LLE and mode coupling 
equations models and Hansson et al. in [[25] showed that both models can be numerically solved 
in similar computational times. 

The generalized mean-field LLE equation is: 
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which is the nonlinear Schrodinger equation (NLSE) with damping, detuning and external 
pumping, which accurately describes Kerr comb generation. In this equation E(t 1 t) is the 









complex envelope of the total intracavity field and hereafter simply the held, t is the so called 
slow time variable, T is the fast time variable, tR is the round trip time, a is half of the total loss 
per round trip which includes the internal linear absorption and coupling loss. So is the phase 
detuning, L is the cavity length, /j/ is the border dispersion coefficient, y is the Kerr coefficient, 
6 is the coupling coefficient between the waveguide and the microresonator and E in is the pump 
held (normalized such that P m = C m | 2 ). To facilitate our analyses in a more general framework 
we will additionally use the normalized detuning and pump held according to Eq. 2 and 3. 
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Recently, some authors explored and identihed regions corresponding to different types 
of operation in the detuning and pump power (A, l^l 2 ) parameter space in the anomalous 
dispersion regime |I 6 ] 19 23). The characterized regions are: stable modulation instability 
(SMI), unstable modulation instability (UMI), stable cavity solitons (SCS), unstable cavity 


solitons (UCS) and continuous wave (CW). Additionally, Erkintalo and Coen in [221 explored 
the hrst-order coherence properties of each region. 

Single CS generation is a desired goal in much research because it yields a high coherence, 
single free spectral range (FSR) Kerr comb. It has been suggested that this style of comb can 
be used as information carriers in optical communications and optical memories. Although 
the hrst experimental observation of a single CS in microresonators was reported in (26), the 
difficulty to control the number of CSs in the SCS region was also emphasized. In particular, in 
the usual approach in which detuning is swept to a predetermined final value with input power 
held constant, the number of CSs generated is probabilistic. In the current work, we show a 
method not only to obtain a single CS with a high degree of certainty avoiding the chaotic and 
unstable regions such as UMI and UCS, but also to compress it, moving the system through the 
SCS to high power values. 


2. Deterministic Single CS Generation 

The LLE is simulated numerically using the split step Fourier method (SSFM). The simulation 
parameters used in this work are: tR = 1 /226 GHz, /T = —4.7 x 10 26 s 2 m *, a = 0.00161 


y = 1.09 W 'm *, L = 2n x 100 jUm and 9 = 0.00064 27[. These parameters correspond to a 


SUN 4 microring resonator of 100 ll m radius with anomalous dispersion, loaded quality factor 
Q of 1.67 x 10 6 and photon lifetime f ph = tn/(2a) = 1.37 ns. Addionally, the simulations are 
done initializing the intracavity field in the frequency domain E(co) to a circularly symmetric 
complex Gaussian noise field with a standard deviation cr noise = 10 2 [W 1 / 2 ], which with 
normalization P = \E\ 2 is equivalent to a mean power of -150 dBm per cavity mode and uniform 
random phase. 

To show the probabilistic nature of the number of CS after a detuning swept at constant pump 
power we perform simulations according to the insets of Fig. 1(a). These insets depict how the 
detuning is swept linearly from A = 0 to a final value of A f over a time interval of 0.3 jis («218 
f ph ) at a constant |S | 2 = 18.9 (P in = 180 mW). Fig. 1(a) shows the number of temporal peaks 
(in the SCS region this corresponds to number of CSs) versus final value of detuning A f . The 
number of CSs obtained depends very sensitively on A f even though the noise initialization was 
the same for all the simulations of Fig. 1(a). Additionally, the number of CSs is different when 
we repeat the simulation with the same detuning and pump parameters but different realizations 
of the initial noise, as shown in the histogram of Fig. 1(b). The behavior of the total intracavity 
energy 6 / trllrLI , which is proportional to the integral of the intensity over a round trip time, give 









us an idea of the Kerr comb stability. Fig. 1(c) shows total intracavity energy during a detuning 
swept at constant input power. The noisy behavior in the yellow segment gives us evidence 
that the system is in the UMI region. Then, the oscillatory behavior in the blue segment is a 
clear characteristic of UCS region. Finally, the total intracavity energy becomes stable when the 
microresonator approaches the SCS region. Our hypothesis in this work is that the uncertainty 
in the number of CSs in the SCS region results from chaotic field evolution while the system is 
crossing the UMI and UCS regions. 





Fig. 1. (a) The number of peaks as a function of final value of detuning Af with A swept as 
shown in the insets. The pump power |S| 2 is set to a constant value of 18.9 ( P m =180mW). 
These simulations were done using the same realization of initial noise, (b) Histogram of 
number of peaks for 1000 simulations with different realizations of initial noise, using the 
same detuning and pump parameters in the simulation of (a) with A f = 12.42. (c) Total 
intracavity energy (blue) and detuning (green) as a function of time for the simulation of 
(a) with A f = 12.42. 


Although some of the boundaries between the regions can be solved analytically (e.g. 
Hamiltonian-Hopf limit between CW and SMI regions (23J), there is still no theoretical 
solution for the boundary between the SMI and UMI regions. Therefore, we characterize the 
(A, |5j 2 ) parameter space via numerical simulations to gain knowledge about these boundaries. 
We initialize the microresonator pump at a detuning A = 0 and power |S| 2 = 6.3 (60 mW) 
with well-defined behavior in the SMI region. This point corresponds to the location of the 
red triangle in Fig. 2(a) where a Turing rolls pattern 121 24j of 13 equally spaced peaks 
is formed. The intensity and power spectrum obtained at this point are shown in Fig. 2(c) 
(left) and (center), respectively. This pattern corresponds to a Kerr comb with frequency 
spacing between lines equal to 13 times the microresonator FSR and possesses low intensity 
noise and high stability, as shown by its steady state intracavity energy in Fig 2(c) (right). 
According to Erkintalo and Coen in |22) , the UMI region exhibits high first-order coherence. 
Coillet and Chembo in |[24j demonstrated that Turing rolls patterns achieve phase locked states 
independently of the initial conditions and can be easily generated by decreasing the optical 
frequency of the pump laser towards the resonance, as presented in their experimental examples. 
These reasons make Turing rolls pattern a good choice for our initial point. After 200 ns, we 
jump in a single step to a point in (A, |S| 2 ) parameter space and keep the system there for 1.8 /is 
(«1310 f ph ). Simulations were repeated using the same realization of initial noise for different 
final values of detuning A and pump power |5j 2 up to 6.21 and 10.5 (100 mW), respectively. 
Fig. 2(a) shows the number of time domain peaks for each final point. Each point (A, |S| 2 ) 
























is classified based on intensity, spectrum and total intracavity energy stability at the end of 
the simulation in order to determine the boundaries of the different operating regions. The 
characteristics of each region are shown in Fig. 2 (c) to (f) and the result of this classification 
is shown in Fig. 2(b). Notably, U Mm in the UMI and UCS regions is characterized by noisy and 
oscillatory behavior, respectively, while in the SMI and SCS regions it is constant as mentioned 
before. 
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Fig. 2. (a) The number of temporal peaks as a function of the final point in (A, |S| 2 ) 
parameter space at low pump power. The red triangle shows the initial point of the 
simulations. These simulations were done using the same realization of initial noise. The 
black curves correspond to CATs using Eq. 4 with |5'|offset °f -0.42, 0, and 0.42 equivalent 
to offsets of -4, 0 and 4 rnW respectively, (b) Characterized regions according to the 
features of each region shown in (c) to (f). (c-f) Final intensity (left), spectrum (center) 
and intracavity energy versus time t (right) for the point (c) (Af = 0, |Sfj 2 = 6.3) in the SMI 
region, (d) (A f = 3.9,|5 f | 2 = 10.5) in the UMI region, (e) (A t = 4.8,|S f | 2 = 10.5) in the 
UCS region, and (f) (A f = 4.8, |S f | 2 = 6.3) in the SCS region. 


Because instantaneously jumping to a single final state in detuning-power parameter space 
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as in the simulations of Figs. 2(a) is not physically realizable, we next define a path that can 
be followed smoothly between our starting point (A = 0, |.S| 2 = 6.3) and a target point chosen 
to fall within the single CS region (A = 4.35, |S| 2 = 5.04). In particular, we select a path that 
goes around the chaotic and unstable regions. We call such a path a chaos-avoiding trajectory 
(CAT). In this study we consider a CAT given by the following: 

|S| 2 = 4.15 exp (—3.09 A) + 2.15 exp(0.196A) + |S| 2 ffset (4) 

This CAT is defined using a curve-fitting tool and takes the form of a two-term exponential 
function. Additionally, we introduce a power offset parameter 5'| 2 )ffs( . t which allows us to 
describe a family of CATs, portrayed by the black curves in Fig. 2(a). To test the performance 
of the CAT we run simulations in three stages. In the first stage we dwell at the selected initial 
point (A = 0, |*S| 2 = 6.3) for 1.0 /is («728 f ph ) to produce a steady state Kerr comb in the 
SMI regime. Next, we sweep the detuning linearly from A = 0 to a final value A f over a 0.3 
/is interval («218 f ph ) while varying the input power according to the CAT in Eq. 4. We then 
continue the simulation at fixed detuning and input power for an additional 2 /is («1445 f ph ). 



Fig. 3. (a) The number of peaks as a function of final value of detuning Af with A swept 
as shown in the inset with pump power adjusted through the CAT, Eq. 4, with |S|Q ffset = 0 
using the same realization of initial noise for all simulations, (b) Histogram of number of 
peaks for 1000 simulations with different realizations of initial noise with the same pump 
power and detuning parameters of the simulation of (a) with A f = 7.45. (c-d) Region in 
which a single CS is generated when a uniform offset |S|Q ffset is applied to the CAT for 
different values of final detuning Af with (c) constant detuning interval of 0.3 /is and (d) 
constant detuning speed of 25 units of A per /is. 


Figure 3(a) shows the number of temporal peaks versus the final value A f for fixed realization 
of initial noise. The insets illustrate the coordinated variation of A and |.S| 2 in time according to 
the CAT. The green area depicts a range of A f values (from 4 to 10.8) for which a single CS is 
always obtained. In another case we keep A r fixed but repeat the simulation for 1000 different 
realizations of initial noise. As shown in the histogram of Fig. 3(b), a single CS is obtained every 
time. Fig. 3(c) is similar to Fig. 3(a), in that we keep the noise initialization fixed and vary A f , 
except that now we use different values of |A 1 1offset from -0.52 to 0.52 equivalent to an offset 
of -5 to 5 mW to the CAT. The green shaded area again shows the cases where we generate a 
single CS. Additionally, we test the CAT performance keeping the sweep rate of the detuning 
constant at 25 units of A per /is for different values of A f . This is in contrast to the simulation of 
Fig. 3(c) in which the detuning in swept over a constant interval of 0.3 /is, implying different 





























sweep rates for different final detunings. The simulation results are displayed in Fig. 3(d), 
which again shows a green area in which a single CS state is reached. Moreover, we perform 
simulations using half and twice the sweep rate of the detuning with very similar results of Fig. 
3(d). All these simulations demonstrate that the CAT for repeatable generation of a single CS 
is nonunique and can be robust over a finite range of operation. It should be emphasized that 
this CAT could take other forms different than a two-term exponential function, provided that 
it avoids the unstable and chaotic regions. 

We note that thermal nonlinearities play an important role in practical microresonators. 
The response time of such thermal nonlinearities is much slower than that of the Kerr 
nonlinearity, which complicates computational studies. For this reason thermal nonlinearity is 
not included in the simulations in this paper. Further research is required to determine whether 
the current scheme for deterministic generation of single cavity solitons remains effective with 
the inclusion of thermal nonlinearity. At a minimum we anticipate that it may be necessary to 
slow down the CAT sweep speed to a time scale consistent with the thermal nonlinearity. 

3. CS compression 



Fig. 4. The number of temporal peaks as a function of the final point in (A, |S| 2 ) parameter 
space at high pump power. The red triangle shows the initial point of the simulations. 

These simulations were done using the same realization of initial noise. The black curve 
corresponds to a compression function using Eq. 5. The initial point of this compression 
function is the final point of the CAT with |S|o ffsct = 0 in Fig. 2(a). (b) Characterized regions 
using the same criteria of Fig. 2(b). 

Our final goal is to demonstrate compression behavior when we move the system from 
relatively low to relatively high detuning and power under a trajectory that remains within the 
SCS region. A process similar to the characterization of Fig. 2 is performed to explore behavior 
over a larger range of parameter space (final values of detuning A and input power |5j 2 up to 
31 and 105 (1 W), respectively), starting at initial point of detuning A = 0 and pump power 
|5j 2 = 63 (0.6 W). This point has higher power than the initial point used in Fig. 2(a) because in 
this higher range of power, the microresonator needs more initial energy to show all the possible 
regions of behavior without collapsing to the CW state. The results of the simulations in this 
range and the classification are shown in Figs. 4(a) and 4(b), respectively. Then, we define a 
compression function from (A = 4.35, |S| 2 = 5.04), the end of a previous CAT trajectory, to 
(A = 27.45, |S| 2 = 105) using a trajectory given by 

|S| 2 = 2.846exp(0.1316A) 


(5) 

















Again, this function is not unique, and a wide range of end points A and .S’| 2 lead to similar 
results, provided that the trajectory remains within the SCS region. 

To provide a complete picture of the generation and compression processes of a single CS, 
we run a simulation in five stages, as shown in Fig. 5(a), which plots the coordinated change 
in detuning and input power versus slow time. Figure 5(e) shows the intensity and spectrum of 
the steady state initial point (Turing rolls pattern of 13 temporal peaks). In Fig. 5(b) the Turing 
rolls pattern fall into a single CS around 1.4 /is. Fig 5(f) shows the intensity and spectrum of the 
generated single CS. In Fig. 5(c) we can observe how the Kerr comb is broadening while the 
system is going through the compression function between 1.6 /is and 2.6 /is. Figures 5(f) and 
(g) show that the CS is compressed from 73 fs to 29 fs full width at half maximum (FWHM) 
pulse width and the bandwidth Aco/2n goes from 3.98 to 10.2 THz. Figure 5(d) shows f/ Intra 
along a CAT leading to single CS generation and subsequent compression. Neither noisy nor 
oscillatory behavior is seen, providing evidence that the system remains in SMI or SCS regions 
for all simulation time and should consequently be highly stable with low intensity noise. 

Although we are not taking into account high order dispersion terms and other nonlinear 
effects (e.g. Raman) that will be present in real systems, we believe these results make clear that 
once a single CS state is reached compression should be possible while maintaining stability 
through appropriate coordinated increases in A and .S' 2 . 
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Fig. 5. Generation and compression processes of a single CS through the CAT and 
compression functions in Eq. 4 and Eq. 5, respectively. From top to bottom: (a) pump power 
(blue) and detuning (green), (b) temporal intensity, (c) spectrum, (d) total intracavity energy 
vs. slow time t. (e-g) Intensity (right) and spectrum (left) at (e) t = 1 /is, steady state initial 
point, (f) t = 1.6 /is, after single CS generation and (g) t = 3.6 /is, after compression. 


4. Conclusions 

In this work, we have presented a novel method to generate a single CS in anomalous 
dispersion microresonators in a highly deterministic way through coordinated tuning of pump 
























































frequency and power in order to avoid chaotic and unstable operating regimes. The simulation 
results presented here could help accelerate progress towards establishing microring resonators 
as highly coherent and stable single FSR Kerr comb sources for practical applications. 
Furthermore, we have also demonstrated a way to compress the single CS by further 
coordinated variation of the pump parameters in the region where CSs are stable. 
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